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Abstract 

Using detailed balance and scaling properties of integrals that appear in the Coulomb gas refor- 
mulation of quantum impurity problems, we establish exact relations between the nonequilibrium 
quantum decay rates of the boundary sine-Gordon and the anisotropic Kondo model at zero tem- 
perature. Combining these results with findings from the thermodynamic Bethe ansatz, we derive 
exact closed form expressions for the quantum decay rate of the dissipative two-state system in 
the scaling limit. These expressions illustrate how the crossover from weak to strong tunneling 
takes place. We trace out the regimes in which the usually applied Golden Rule (nonadiabatic) 
rate expression fails. Using a conjectured correspondence between the relaxation and dephasing 
rate, we obtain the exact lower bound of the dephasing rate as a function of bias and dissipation 
strength. 
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Quantum impurity problems (QIPs) have attracted a great deal of interest recently. This 
is because the underlying physics is non-trivial and the models are manageable technically 
despite their essentially nonperturbative nature. In addition, they have a multitude of 
experimental applications, including the Kondo effect, quantum dots, dissipative quantum 
mechanics, tunneling in quantum wires and fractional quantum Hall devices ||T|. 

There have been discovered various relations between thermodynamic quantities of the 
anisotropic Kondo model (AKM) and the boundary sine-Gordon (BSG) field theory model 
0. Each of these integrable models is of considerable interest and it is remarkable that they 
are closely related. The AKM corresponds to the dissipative two-state system (TSS) in the 
Ohmic scaling limit and the BSG model corresponds to the Schmid model ||]. The latter 
describes a quantum Brownian particle coupled to an Ohmic heat bath and moving in a 
tilted cosine potential. The tight-binding (TB) limit of the Schmid model is equivalent to 
the strong-backscattering limit of the BSG model, while the weak-binding (WB) limit of the 
Schmid model represents the weak-backscattering limit of the BSG model. In the sequel, 
we use the language of dissipative quantum mechanics. 

The equivalence or difference of these models can most easily be seen in the Anderson- 
Yuval Coulomb gas representation for the partition function. The partition function of all 
these models can be expressed as the partition function of a one- dimensional classical gas of 
positive and negative unit charges with "log-sine" interactions. The AKM (or TSS) problem 
differs from the BSG (or Schmid) model by a different ordering prescription. 

In this Letter, we calculate the full nonequilibrium quantum decay rate of the TSS at 
zero temperature in closed analytic form. Upon using a conjectured correspondence between 
relaxation and decoherence in the TSS, we also determine the lower bound for decoherence 
in this system. These outstanding problems are solved by establishing exact functional 
relations between nonequilibrium quantum rates of the TSS and the Schmid model. 

The TSS Hamiltonian in pseudospin form reads 

H TSS = -lA T <j x -l[e + f{t)}a z , 

where h = k B = 1. Here A T is the tunneling coupling, e is the bias energy, and f(t) is a 
random force with Gaussian statistics. All effects of f(t) are captured by the second integral 
of the force autocorrelation function, denoted by Q(t). In the scaling (or field theory) limit, 
we have Q{t) = 2K ln[ (oj c j nT) sinh(7rT|t|) ] + i nK sgn(t), where T is temperature and u c is 
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a cut-off frequency. The Kondo parameter K is a dimensionless Ohmic damping strength, 
which is the inverse of the Luttinger parameter g or the filling fraction v in quantum impurity 
problems. 

The TB Schmid model in second quantized form with coupling energy A s and bias energy 
e is given by 

H s = -| A sXl ( a n a n+i + h -c-) - [e + f(t)]J2nala n . 

n n 

The nonequilibrium dynamics of both models may be computed using the Keldysh or 
Feynman- Vernon formalism for the reduced density matrix (RDM). Consider a particular 
path on the (q, g')-plane of the RDM parametrized by charges {uj = ±1} and {vi = ±1} 
chronologically arranged on the q and q' path, respectively. The influence functional due to 
the stochastic force f(t) for k charges on the g-path and i charges on the g'-path, divided 
up into the self-interactions of the paths and the interaction between the paths, reads 

T = exp < UjUiQ(tj - U) 
< j>i=i 

i e k n 

+ £ vjvtQ*^ - o - j2 E -**)}■ 

j>i=l j=l i=l * 

For the TSS the charges Uj,Vi alternate in sign, while for the Schmid model they may be 
ordered arbitrarily. At long times, the Schmid model is a Poissonian transport model: the 
population dynamics results from uncorrelated direct forward and backward transitions by 
n wells, n = 1, 2, ■ • ■ . Since they are not classical, the respective weights per unit time 
("rates") 7^ are not necessarily positive because of quantum interference 0. In the TB 
representation, the Fourier transform of the population distribution P(k,t) = ^ n e lkn P n (t) 
takes the form P(k, t) = Yl™=i ex P [t { e%kn ~ l) In + 1 (e~ lkn — l) 7"] , where the rates 7^ 
describe uncorrelated direct forward/backward transitions by n TB states. 

The Poissonian dynamics is found upon carrying out a cluster decomposition of the 
individual path contributions to the Laplace transform P(k, A) in the limit A — > 0. A cluster 
is a A-independent (irreducible) path section which starts and ends in a diagonal state of the 
RDM, Ui = Ylj v j an d cannot be factorized into clusters of lower order. By definition, 
the clusters are noninteracting. Therefore a succession of clusters factorizes and the integral 
over a (time) interval separating any two clusters produces a factor 1/A. Path segments 
which start and end in a diagonal state and again have interim visits of diagonal states 
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are reducible, i.e. they factorize into clusters of lower order and factors of 1/A. However, 
after subtraction of all the reducible components, an irreducible part is left. The irreducible 



the transition rates 7*. 

The irreducible parts of arrangements with 2(n + i + k) charges contribute to the order 
£2(n+t+k) Q f They can divided into a multitude of different subsets. Each subset a 
is a particular arrangement of n + 21 time-ordered u-charges and n + 2k time-ordered v- 
charges. The number of subsets with different order of the u- and f-charges is N = N U N V = 



arrangements. Each subset a defines a partial forward rate 7"' + . The corresponding partial 
backward rate 7™'" is represented by the arrangement with reverse order and reverse sign 
of the charges. The individual partial rates are directly related by detailed balance, 7"'" = 
e~ ne l T 7^' + . Detailed balance is satisfied for every partial rate because the interaction Q(z) 
between the u- and f-charges is analytic for complex time z in the strip > Imz > — 1/T 
and has the property Q(t — i/T) = Q* (t) . 

In the remainder, we restrict the attention to zero temperature. At T = 0, all frequencies 
of the Schmid and TSS model can be absorbed into a single dimensionless parameter, x s = 
{e/uo c ) K A s /eandx T = (e/u c ) K A T /e. For later convenience we also introduce the frequency 
scale e (analogous to the Kondo scale in the Kondo model and the scale Tg in QIPs [0) 
and the dimensionless bias v = e/e . The scale e is related to the bare parameters of the 
Schmid and TSS model as 

" T 2 (K) uo 2 c K ' 2 2K u 2K ' 

In the scale e , the selfduality of the BSG model becomes evident. In the second form, we 
have anticipated from below, that at fixed e the TSS bare coupling A T is related to the 
BSG bare coupling A s via A T = 2 sin^if) A s 0. 
The rates of the Schmid model can be written as || 



To proceed, we map the u- and f-charges onto a single time axis. To this end, we introduce 
charges rjj = ±1 describing forward/backward moves along the quasiclassical path q + q' 
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contributions of all charge arrangements with excess charge ^ ■ Uj = ±n, ^ Vi = ±n define 





(1) 



£=n 



and charges £j = ±1 representing moves along the quantum fluctuation path q — q'. The 
cumulative charge pe = Yll=i £fc measures how far the system is off-diagonal after £ moves. 
Then we have 

Ki = E E [ cos ^ # ± sin ^ 7 *1 • ( 2 ) 

Here {•■■}' denotes the constraints = for the 2£ ^-charges and YliVi = 2n for 

the 2£ 77-charges. The number of different charge sequences contributing to U^ e is N^ v = 
(e+n)f(e-n)\ • ^ ne i ma g mar Y parts of the bath correlations in the influence functional T 
add up to the influence phase £ (£,f/) = i\K Y^=i PjWj- 

All the quantum fluctuations are in the 2£— 1-fold integrals if, which are independent of ff. 
With the notation J °° T>it-\(j) ■ ■ ■ = Jo^ ^ r 2£-i ' ' ' ■ ■ ■■> where the Tj are the dimensionless 
charge intervals (scaled with e), we get 



PCX) 

It (15 ± i II (15 = / ^-i (f) Gf (f , a 

./0 



The bias phases are combined in the phase (p e (r,£) = Y^Li Pj T j- The interactions of the 



^-charges are encapsulated in the factor (we put r-j = tj — U — J2k=i T k ) 



2£ 



Gf\r,0 = exp (2KJ2 fcfchr,,) - ^(^0 • 
j>»=i 

The subtractions are to ensure that the if represent irreducible clusters. Subtractions are 
needed whenever one or several of the cumulative p-charges are zero. For instance, putting 
£ = 2 and & = -fi, U = we then have Gf ubtr) = (7^3) - 2 *. 

The expression for can be simplified considerably upon utilizing the following obser- 
vations: 

(i) At T = 0, backward moves to higher wells are absent, i.e. the aforedescribed partial 
backward rates 7"'" are all zero. Hence for any subset a, the respective linear combination 
of the //-integrals cancels that of the /^"-integrals, while they add up in 7"'" 1 ". 

— # 

(ii) For every individual ^-configuration, in which all cumulative charges p- have the 

I — # 

same sign (some of them may be zero), the integrals l t are simply related by = 
tan(27rfT£)/7(0. 

Property (i) is of general nature since it is due to detailed balance. Property (ii) is 
based on the scale invariant logarithmic charge interaction in and is responsible for 



5 



integrability. It does not hold when scale invariance is broken by finite cutoff or finite 
temperature. 

Upon exploitation of properties (i) and (ii) a multitude of linear functional relations 
between the various integrals with different £ but same i can be derived. The analysis 

yields that all integrals Ig(£) with one or several equal to zero can be expressed in 
terms of those which have all p nonzero. Making use of these relations we find that there 
are formidable cancellations among the various contributions to U^ t In particular, the 
contributions from arrangements with any fixed nonzero number of negative ^-charges cancel 
each other exactly. Thus we have U^g = for i > n. Hence the contributions of all paths 
which undertake backward moves between virtual intermediate states add up to zero and 
only the paths with the minimal number of flips 2n contribute to the forward rate 7+. We 
finally get 

Kn = E J n (3 II M«KPk) ■ (3) 

We now turn to the study of the incoherent dynamics in the TSS. The subset of charge 
sequences in the Schmid model in which both the u- and f-charges alternate in sign defines 
a partial rate of r y 1 . This partial rate, denoted by t^, is the full rate describing incoherent 
relaxation in the biased TSS. At T = 0, the expansion in the number of transitions (fugacity 
expansion) gives 

00 

T ± = ^2 Iff ; 7? = «#W?, (4) 
1=1 

where W i = U 1 g , and where the double sum in @ is subject to the TSS constraints. The 
77-sum is easily performed, yielding 

I _ COS \7rK)J2 Vt ± Ci tan(vrir)/-] . 
V 7 {^}" 

Here {•}" denotes the TSS constraint £,2k = —(,2k-i- 

Making again extensive use of properties (i) and (ii), we find that all the partial backward 
rates 7" vanish and that can be expressed in terms of the integrals, in which all 

Pj are nonzero. In the end we find using (^) 

ff M- i r I ^ff;.. (5) 

[2 sm(7rA )\ zn 



w~- 
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Employing relation (|5]), we can express the partial rate 7+ of the TSS in terms of the rate 
7+ of the Schmid model. At fixed renormalized coupling e we have 

7+ = (-ir 1 4sin 2 (^)7+. (6) 

Upon combining the thermodynamic Bethe ansatz with the Keldysh approach, the rate 7+ 
has been found as 

7 + = LIT! r(j)r(*rn) j_ {2K _ 2)n 

Jn n\ r[| + (K- l)n] 2tt 1 J 

Use of the relation @ and the expression (0) gives the partial rate 7+ in explicit closed 

form. The perturbative expansion of the full rate 7 + of the TSS is found to read 

00 

7 + = ^E fl m W^ 2)m - (8) 

* m=l 

1 r(ATm) [1 - cos^vrifm) ] 



a m( K ) 



ml r[|+(A--l) 



m 



Expression (H) is the weak-tunneling expansion of the TSS. For rational values of K, the 
series (§) is a linear combination of hypergeometric functions. 

In the regime K < 1, the perturbative series is absolutely converging for large enough v, 
i.e. large enough bias. For K > 1, (^) is absolutely converging for small enough v, i.e. small 
enough bias. The leading term is the Golden Rule rate 7q R = [ii/2T(2K)]ex^ |]. ||. 

An integral representation for 7 + is found by writing T{Km)/T[% + (K — l)m] as a 
contour integral. Within the radius of convergence of series (§) we can interchange the order 
of summation and integration, yielding 

~ 1+ = Re ^nJ c T [{Z ~ 1 ~ ~ {Z ~ 1 " zKUl)h] ' 

where u 1 = v 2K ~ 2 and u 2 = e l27rK u 1 . The contour C starts at the origin, circles anti-clockwise 
around the square-root branch point which is near 1 + u l j 2 for small v 2K ~ 2 , and returns to 
the origin. The integral converges for all K and v, not just in the regime where the series 
(§) does. Therefore it must also yield the asymptotic expansions for small (large) enough v 
when K < 1 (K > 1). 

The asymptotic expansion for K < 1 can be found by changing variable to y — 



[u 1/2 z 



)V(i K ) a nd expanding the integrand for small v. The strong-tunneling series is 

00 

* + =2^£ W** 1 " 1 - ( 9 ) 

* n=0 




FIG. 1: The scaled rate 7 + /7gr is plotted versus x T for various K < 1. The circles represent 
the weak-tunneling series and the squares the strong-tunneling expansion. The full curve is the 
hypergeometric function expression 



For | < K < 1, the branch cut of the M 2 -term is not encircled by the contour. In this regime 
we obtain 

K(K) = d n (K) = i- fl [[ " h -r i = • (10) 



n! (I 



l-A'J 



For K < i both the u x - and « 2 -term contribute, giving 



b n (K) = 2sm 2 [^(l-n)]d n (K). 



(11) 



In the limit K <C 1, the perturbative series (§) as well as the asymptotic series (|9|) can 
be summed to the known result 7 + = ttKA^ / \/ A| + 



1 . 



The asymptotic series (|9]) is an expansion in powers of e 2 / e o- The leading term is inde- 
pendent of e and represents the forward rate for the symmetric TSS, 



! 



(12) 



where b Q (K) is given in fllCf ) and (jTT|), respectively. 

In Fig. |I|, the scaled rate k + = 7 + /7gr is plotted for different K < 1. The horizontal line 
is the K = 1 result. For K < |, the various tunneling contributions interfere always destruc- 
tively, which leads to a reduction of the rate. For | < K < 1, the tunneling contributions 
interfere constructively for all a; T . In the range \ < K < ~, the rate goes through a maxi- 
mum. This reflects constructive interference at small and intermediate x T , and destructive 
interference at large x T . 



The nonperturbative result (12) for the relaxation rate may be compared with the deco- 
herence rate 7 dec describing damping of the coherent oscillations in the TSS 0. This rate 



S 




. K = 3/2 



0.75 
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FIG. 2: The same scaled rate, but now for various K > 1. 



has been calculated for the symmetric TSS in Ref. within the framework of integrable 
QFT. In our notation, the result is 7 dec = 7o~/2 for K < | and 7 dec = sin 2 [^f^y] 7 o~ f° r 
| < K < ~. It is known that the same relations hold for a biased TSS in the regimes K <C 1 
and K close to ~ JT|. Upon combining these relations with Eqs. @ - ( ]TT|) it is natural to 
conjecture that the decoherence rate of a biased TSS is 



7de 



^sin 

n=0 



2r Trif (1 
il-K\2 



-n)]d n (K)v* 



2n-l 



(13) 



in agreement with all the known results in special cases. This expression gives the lower 
(T = 0) bound for quantum decoherence in the entire regime < K < |. 

Finally consider the damping regime K > 1 for large v (large bias). By defining the vari- 
able t = e~ t ' K u 1 i 2 z K , we may expand the resulting integrand in powers of Uyl K . Introducing 
the sections K = p + k with p = 1, 2, • • • and < k < 1, we obtain the asymptotic series as 



r 



c m ( K ) 



(2/K-2)m 



V 771=1 

(-l) m 2T(f )sin[±±£m7r] sin(fm7r) 



(14) 



m! Kr[| + (i-l)m] 

This asymptotic series (0) bears a strong resemblance with that of the self-dual Schmid 
model @, H], 0. The powers follow from those of the perturbative expansion ([|) by the 
substitution K — > 1/K, as also do major parts of the coefficient c m (K). However, it differs 
from self-duality by the alternating sign and by the sine-factors. 

Fig. gives plots of the scaled rate k + for K > 1. The rate goes through a maximum 
which is shifted to higher x T with increasing K. At large enough x T , the strong-tunneling 
contributions interfere always destructively. 
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